Coulomb Blockade Peak Spacing Distribution: The Interplay of Temperature and Spin 
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We calculate the Coulomb Blockade peak spacing distribution at finite temperature using the 
recently introduced "universal Hamiltonian" to describe the e-e interactions. We show that the 
temperature effect is important even at fcsT^O.lA (A is the single-particle mean level spacing). 
This sensitivity arises because: (1) exchange reduces the minimum energy of excitation from the 
ground state and (2) the entropic contribution depends on the change of the spin of the quantum 
dot. Including the leading corrections to the universal Hamiltonian yields results in quantitative 
agreement with the experiments. Surprisingly, temperature appears to be the most important effect. 



PACS: 73.23.Hk, 73.40.Gk, 73.63.Kv 

Among the unique features of quantum dots (QDs) 
the possibility to control their number of electrons N.\ 
This is done by weakly coupling a QD to source and 
drain leads and using a gate voltage Vg to control its 
electrostatic potential. When the thermal energy ksT is 
smaller than the charging energy Ec required to add an 
electron to the QD, the electron transport is blocked and 
N is fixed. By sweeping Vg this Coulomb Blockade (CB) 
effect can be overcome at the particular value where 
the transition occurs. The conductance G{Vg) 

shows then a series of sharp peaks as a function of Vg . 

At sufficiently low temperature, fcsT <C S where S is 
the energy gap between the ground state (GS) and the 
first excited state of the QD, only the former contributes 
significantly to the conductance peak. In that case, the 
position of the CB peak is proportional to the change 
of the GpS energy of the QD upon the addition of one 
electron.ld Therefore, the CB peak spacing distribution 
(PSD) yields information about the many-body GS prop- 
erties of the QD. P 

This has. Jbeen the subject of experimentalETB and 
theoreticalQiij work over the last years. One reason is 
that the simplest model used for the CB conductance 
peaks fails drastically in describing the observed PSD. It 
assumes a constant e-e interaction (Ec) — hence the name 
constant interaction (CI) model. Recently, however, it 
has become clear that residual e-e interactions (i.e., those 
beyond Ec) platan important role in determining the 
GS of the QDCrtJ and therefore must be included in the 
description of the PSD. In papticuiar, the addition of the 
average exchange interactionJ13~b3 which gives the "uni- 
versal Hamiltonian", hereafter called the constant ex- 
change and interaction (CEI) model, leads to a completely 
different PSD.EJ Yet, this is not enough to account for the 
observed distribution: other (smaller) contributions-=j 
such as the "scrambling" effect ,H "gate" |perturbationsE£l 
and the fluctuation of the interactionsEJ — have to be 
consideredJlj So far, however, a much simpler effect has 
not been considered within the CEI model: the effect of 
finite temperature. 

The goal of this work is twofold. First, we show that in 
the CEI model the temperature effects are more impor- 
tant than in the CI model and that they become signif- 



icant even at ksT ^ O.lA where A is the single-particle 
mean level spacing. In particular, the shape of the peak 
spacing distribution changes significantly while increas- 
ing temperature. Since most experiments were done in 
the regime /csT^ 0.3-0.5A — an exception is Ref. || — our 
results are crucial for interpreting the experimental data. 
Second, we calculate the PSD including all the leading 
order corrections to the CEI model. The final result for 
the distribution is in quantitative agreement with the ex- 
perimental data of Refs. ^ and |l^ once the temperature 
effect is included. Surprisingly, the latter introduces the 
biggest correction to the T = CEI model result. 

Mesoscopic fluctuations associated with single-particle 
properties of chaotic QDs are known to be well described 
by random matrix theory (RMT) in an energy window 
up to the Thouless energy i?th- The treatment of the e-e 
interaction is more subtle however. From Fermi-liquid 
theory, we expect the screening of the Coulomb interac- 
tion to be important for TV ^ 1. In that case, the residual 
interaction should be weak, and perturbative treatments, 
such as RPA, seem adequate. Using such an approach 
and RMT to describe .the single-particle Hamiltonian, it 
is possible to deriveE3E3 an effective Hamiltonian Hqb 
for the QD. The small parameter in the perturbation the- 
ory is oc N~2 ^ 1 with g = Eth/A the dimensionless 
conductance. The zeroth-order term {g—>-oo) in this ex- 
pansion co rresp onds to the "universal Hamiltonian" and 
is given bytHltJ 
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where {£„} are the single-electron energies, JV = CgVg/e 
describes the capacitive coupling to the control gate, S is 
the total spin operator, and Js is the exchange constant. 
The difference between the CEI and CI models is the ad- 
ditional term proportional to S'^. Since J5 < A/2 has a 
fixed value, the mesoscopic fluctuations in the spectrum 
of i?QD arise only from This is a key point for un- 

derstanding its GS: while in the CI model the levels are 
filled in an "up-down" scheme — which leads to a bimodal 
PSD — in the CEI model it is energetically favorable to 
promote an electron to a higher level and gain exchange 
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FIG. 1. CB peak spacing distribution obtained using the 
CEI model (Eq. |l|) with T = and Js = 0.32A (r^-LS). The 
dashed (dotted) line shows the A'^ even (odd) contribution 
(each spacing corresponds to a A'^— 1^ A'^— >A'^-|-1 transition). 
Notice that the strong even/odd effect predicted by the CI 
model is significantly reduced and that the 5-function in the 
odd distribution persists — but is shifted by Js- The zero in 
the horizontal axis corresponds to 2Ec in the CI model. 



energy whenever the spacing between two consecutive 
single-particle levels is smaller than 2Js (for N even). 
This leads to a GS with -S'^-L aod the simple "up-down" 
filling scheme bpaks down.ll3ilj Consequently, the PSD 
is very differentEJ from the CI model result (see Fig. 1). 

The T = distribution shown in Fig. 1 will remain a 
good description so long as the contribution from the ex- 
cited states can be ignored. To estimate that, let us cal- 
culate the average occupation V of the first excited state 
assuming that only it and the GS are relevant. In the CI 
model, the energy gap S between those states is (5 = As 
where Ae is the single-particle energy spacing between 
the top levels. Using the CUE Wigner-Dyson distribu- 
tion for Ae — broken time-reversal symmetry is assumed 
throughout— we find V = 0.02 for fc^T = 0.1 A. Thus, 
the excited states can indeed be ignored for fc^T^O.lA. 
This is not the case in the CEI model, where S gets re- 
duced by the exchange interaction. For instance, for N 
even, ,5=|A£-2Js| and "P^ 0.14 with Js = 0.32A. There- 
fore, exchange not only modifies the T = distribution 
but makes the temperature effect stronger. 

Furthermore, at finite temperature, the peak position 
involves the change in free energy of the QD upon adding 
a particle. Then, the spin degeneracy should plpjj pi im- 
portant role through the entropy contributionJljEil Let 
us consider the regime F <C fcsT, A <C Ec, where F is 
the total with of a level in the QD. Near the CB peak 
corresponding to the N-Yr^ ^ ^ transition, the linear 
conductance is given bycJ'E3 
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with Ta the partial width of the single-particle level 



a due to tunneling to the left (right) lead and 

„ = ^Feq(j|iV)|(vI/f|ct_JvI;f-i)f [l-/(,^._e,)]. (3) 



Here Pj^ is the equilibrium probability that the QD con- 
tains N electrons, iJQD|*f) =ej|*f), i^eqOW) is the 
conditional probability that the eigenstate j is occupied 
given that the QD contains N electrons, and /(e) — 
{1+ exp[(e — Since near the peak only the 
states with A^ — 1 and A^ electrons are relevant, we have 
PjY ~ f{TM_~J-N-i) with Tn the canonical free energy 



of the QD.B To make the dependence on M explicit, let 
us denote by {E^} the eigenenergies of ^^qd without the 

charging energy term. Then, €j—ei = E^ —E^^^+2Ec5M 
and Tn-Tn-i = ^^N+'iEc5N with 5N ^ (iV-i)-7V and 



^JiN^Ef-E^'^ + kBT\n[F,^{:j\N)/F,^{^\N-l)] for any 
i and j. The contribution of the transition i— >j to the 
conductance reaches its maximum when 
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In the particular case where the transition between GS 
dominates, and taking the spin degeneracy into account, 
the CB peak position is given by 
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We see that the peak is shifted with respect to its po- 
sition at T = by an amount depending on the change 
of the spin of the QDJl5i22l Because the r.m.s. of the 
PSD is~0.3A (see Fig. 1), this shift is significant even 
for fcsT^O.lA. In addition, while in the CI model this 
introduces only a constant shift between the even and 
odd distributions, in the CEI model it changes the shape 
of both distributions since different spin transitions con- 
tribute to each one. 

Note that because jafj-this shift, the on-peak conduc- 
tance is renormalized.t^Ej Since different spin transitions 
lead to different renormalizations, the average conduc- 
tance peak depends not only on the average coupling to 
the leads but also on Js and on the statistics of the ^ec- 
trum. This explains the small deviations observecS at 
low temperature from the values predicted in the absence 
of the exchange interaction (in particular the increase in 
Fig. 2 of Ref. H). 

In the general case more than one transition con- 
tributes to the conductance, and the CB peak position 
must be determined by maximizing Eq.(||) with respect 
to Af. For simplicity, we restrict ourselves to the case 
where only the GS and first two excited states are rele- 
vant. By comparing the peak spacing distribution with 
and without including the second excited state, we found 
that this is the case for ksT < 0.2A, which is close to 
the experimental regime. 



2 



a 
Q 

O 



3 - 



1- 



T = 0.1A 



-0.5 




1.5 



G 

Q 



2 



I r.m.s. = 0.24 | 



T = 0.2A 




-0.5 0.5 1 

spacing [A] 



1.5 



FIG. 2. Finite temperature CB peak spacing distribution 
obtained using the CEI model [Eq. @]. Here Js = 0.32A and 
fcsr^O.lA (0.2A) in the top (bottom) plot. Notice that the 
5-function is smeared out and shifted by fcsrin2 and that 
the even distribution develops a peak. 



Figure 2 shows the PSD for non-zero temperature. The 
value Js = 0.32A corresponds to a gas parameter r^. ~ 1.5. 
As expected, the sharp features are smeared out by tem- 
perature, including the (5-function in the odd distribu- 
tion. But more important, they are shifted due to the 
entropic term in (^). This is particularly clear for the 
peak associated with the ^-function: it occurs at Js in 
Fig. 1 but at Js + kBT\n2 in Fig. 2 (the (5-function is 
always associated with the spin transition — > ^ ^ 0). 
There are two other important effects worth comment: 
a) the even distribution develops a peak at small spac- 
ings while the odd one gets broader — in particular, at 
ksT = 0.2A the maximum of the total distribution is 
dominated by the even distribution, in contrast to what 
occurs at T = 0; b) the relative weight of the long tail of 
the even distribution is strongly reduced and the distri- 
bution becomes more symmetric. 

The peak in the even distribution arises from cases 
where S=l and 5* = states are (almost) degenerate — 
it corresponds to the sharp jump at zero spacing in 
Fig. 1. According to Eq. (|) the CB peak is shifted by 
ifcsrin(4/2) in that case, which gives a shift of A;BTln2 
for the peak in the PSD. The deviation from the T = re- 
sult is still noticeable for fc^T ^ 0.05A (data not shown) 
which is the temperature in Rcf. |[ We found that the 



r.m.s. of the total distribution decreases monotonically 
when increasing T in the range fc^T < 0.2A, in contrast 
to the results obtainecLfor the CI model.Ea The avail- 
able experimental datafl cannot discern this difference. 
An expected effect of the temperature is to increase the 
probability that S ~1. This can be easily estimated by 
calculating the average probability {Fcq{S = 1|A^)). We 
find (Foq(S'=l|7V)) = 0.21(0.31) for kBT = (O.lA). 

So far we have considered only the mean value of the 
residual interactions. It is clear that even at this level 
of approximation, the finite temperature PSD is quite 
different from the widely used T = CI model result: we 
should expect only a weak even/odd effect or asymmetry 
in the experimental data for kBT>0.2A. 

Nevertheless, the distribution still does not agree with 
the one observed experimentally. Thus, we must go to the 
next level of approximation and include the leading order 
corrections to Hqt). There are three contributions that 
lead to corrections of order A/y^ to the spacing: 1) the 
"scrambling" of the spectrum when adding an electron to 
the QD that originates in the change of the charge dis- 
tribution in the QDp 2) the change in the-single-electron 
energies when the gate voltage is swept?l3 3) theJluctu- 
ation of the diagonal part of the e-e interaction.EiLAfter 
including them, the Hamiltonian of the QD readsta 
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QD — ^^QD ■ 2 

Q,/3,7,(5 

+ E ^1-^/5,. [{n ^f)xlp+s^^xlp] (6) 

a,/3,cr 

where, because of the fluctuations of the single-electron 
wavef unctions, the matrix elements H^jf}^ g and ^ are 
gaussian random variables. They are characterized by 

var(A'^_^) = 6,, AV5 and var(i?(^/»^ ,) = C2A2 ln(e25)/ff^; 
their mean values can be included in the definition of Ec 
and Js, so that {H^^lp'^g) = {Xlp)^0. Only the diago- 
nal terms of the residual interaction (i.e. those where the 
operators and c are paired) are included in Eq. (^) — 
the off-diagopaLterms can be neglected when calculat- 
ing the PSD.tJlij Here, 5U is taken with re spec t to the 
state with N electrons. We use g = 0.384^iV/2 for the 
dimensioiiless conductance, which corresponds to a disc 
geometry.E3l3 The coefficients hjj , C2 '<md £2 are geometry- 
dependent. We estimated &00 ^ 0.01,Eil C2 -0.003, C2 ~ 37 
for the experiment in Ref. |[ and we assume 611 = 600 • No- 
tice that the magnitude of the "scrambling" and "gate" 
effects [last term in Eq. (||) -axe. much smaller than used 
previously in the literature.c3'L3 

Figure 3 shows the PSD including these corrections to 
the CEI model. The additional fluctuations increase the 
smearing of the remaining pronounced features of the dis- 
tribution as well as of its r.m.s.. The distribution is less 
asymmetric but the even/odd effect is still noticeable — it 
should be kept in mind that the experimental noise will 
contribute significantly to this smearing. At A;bT = 0.2A 
the peak of the distribution is still dominated by the even 
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FIG. 3. Finite temperature CB peak spacing distribution 
obtained from the Hamiltonian (Eq. ^ with A'^^500 {g!v6), 
Js =0.32A and fcsT^O.lA (0.2A) in the top (bottom) plot. 
Notice that both the even/odd effect and the asymmetry are 
small. 



distribution. It is important to emphasize that this par- 
ticular feature is exclusively related to the temperature 
effect. A doteiled analysis of the experimental data of 
Ref. H showal3 that the least noisy data present this sig- 
nature and that the r.m.s. is of order of 0.25A_which 
is consistent with the prediction of this approach.Ea Fur- 
thermore, it is clear from the figures that the temperature 
effect is the main cause of the deviation from the CEI 
model result. 

In conclusion, we have shown that the presence of the 
exchange interaction imposes a more restrictive condi- 
tion on T for observing GS properties in QDs. In par- 
ticular, most of the experiments done so far require in- 
cluding temperature effects for their interpretation. The 
observed PSD seems to be the result of the addition of 
several small contributions. 
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